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1. introduction 

The Liouville and the matrix- mo del [Q] descriptions of 2D quantum gravity^, known 
also as 2D string theory, are complementary to each other in the sense that the each of 
the two approaches has its advantages for certain class of problems. The correspondence 
between these two approaches is historically the first and technically the most elementary 
example of the matrix-string duality. 

In the theory of closed strings (world sheet without a boundary), the predictions of the 
Liouville theory are restricted to the exponents characterizing the scaling behavior of the 
correlation functions. These have been formulated long ago in |^ and and are known as 
the KPZ rule, relating the "flat" and the "gravitational" dimensions of the matter fields. 
The KPZ rule was checked on various microscopic realizations of the 2D quantum gravity 
in terms of statistical models on random lattices and in the last years have been applied 
successfully to evaluate some difficult "flat" critical exponents 0. 

The situation is much more interesting in the case of a world sheet with boundaries, 
which describes open string amplitudes. The simplest such manifold is the disc. In this case 
one can compare not only critical exponents, but also functions of several dimensionless 
quantities. Indeed, an n-point boundary correlator depends on n dimensionless parameters, 
which can be made of the area of the world sheet and the n distances between the points 
along the boundary. 

A well known example is the bulk one-point function in presence of a boundary, 
which plays an important role in 2D quantum gravity [0J^. This function depends on 
the dimensionless ration \ibI ^fj^ and correctly reproduces the results obtained using the 
statistical mechanics of randomly triangulated surfaces. 

Another, much less trivial example is the boundary 2-point function for two boundary 
operators intertwining between the free (Neumann) and fixed (Dirichlet) boundary con- 
ditions. This function depends on two dimensionless parameters associated with the two 
boundary cosmological constants along the Dirichlet and Neumann pieces of the boundary 
and have been first calculated in [§] using the representation of the matter fields as a gas 
of non-intersecting loops (domain walls) on a randomly triangulated surface. 

Other solvable examples of boundary correlation functions are given, again in the 
loop gas formulation, by various configurations of domain walls connecting points at the 
boundary. Such correlators have found application in some statistical problems associated 
with percolating clusters and polymers [If. 

Recently, due to the impressive results in [|10[ and , it became possible to compare 



the results of the matrix model for the boundary correlators with the predictions of the 



Liouville theory. The authors of calculated the general boundary two-point correlation 
function, and the boundary three-point function was calculated in ||1 1|| . 

In this paper we will show that the expression for the boundary Liouville two-point 



correlator obtained in [jT^ is in perfect agreement with the results obtained using the 
microscopic definition of 2D quantum gravity as a loop gas on a random planar graph. 
Moreover, we will give a simple geometrical meaning of the remarkable functional equation 



for the two-point correlator, obtained in |]Tn . 



^ A good review of both approaches is given in Q . 



1 



2. The boundary 2D gravity as a boundary Liouville theory 

2.1. World-sheet description 

The world-sheet formulation of 2D quantum gravity (see, foe example the review 0) 
involves a Liouville field ip and a matter field defined on the world-sheet manifold. From 
the point of view of the 2D string theory these two fields represent the two components of 
the position field of the string. It is most convenient to choose the world sheet with the 
disc topology as the upper half- plane in which case all the curvature is concentrated 
at the infinite point. Then the action is given by the integrali 

Af,x]=j d^z(^£-(y^'^ + Vx^) + t^e^'')+ j dx fiB e"^ + ghosts (2.1) 

and the background charges are defined by the asymptotic of the fields at |x| ^ oo 

(p {g + l)logzz, X (fi'-l)log2^ (2.2) 

The central charges of the matter and Liouville fields are 

c^ = l + b , c^ = l-b (2.3) 

so that + = 26. 

2.2. Bulk and boundary states 



The bulk and boundarya KPZ states are marginal operators Vp and Bp representing 
products of Liouville and matter vertex operators 

Vp(2, z) = e^9+l-\p\M^,z) ^-^{p-9+l)x{z,z) ^ ^2.4) 

Bq{x) = e(4i-l9l)^(^) e-^(9-^)x(x) ^2.5) 

Here we use notations, which are more natural from the point of view of string theory, 
where p make sense of target-space momentum. The dimensions of the matter and the 
Liouville components of the bulk KPZ state Vp are 

A^VA- ''-^-'^^ AnVJ = ZJ!l±<£±i)! (2.6) 



^ Here we adopt the notations commonly used in the Coulomb gas picture, while in the no- 
tations used in ||lO[| the duality transformation takes simpler form. The relation between our 
notations and these of |1CI[ | is given in Appendix A. 

^ The boundary operators in 2D quantum gravity have been first discussed in |12] 



2 



so that its sum is equal to 1. For the boundary KPZ state these two dimensions read 



A.[B,1 ^ A^BJ ^ Zl£!±k±i)!. (2.7) 

In the 1+1 dimensional string theory, the operators Vp and V_p are interpreted as 
left- and right moving closed string states. Similarly, the operators and B_g describe 
left- and right moving modes of open string states with target-space momentum p. These 
states are related to each other by hermitian conjugation. 

There is another, less trivial symmetry, associated with the Liouville direction of the 
target space. This symmetry exchanges the operators (|2.4|) and the operators with the 
same dimension ( |2.6p , but with the "wrong" sign of the Liouville dressing 

Vp(z, z) = e^^'+i+IPl)^ e-*(-P-9+i)x^ (2.8) 

B,(x)=e(^+l^l)^e-^(-'^-"^^)^. (2.9) 

The correlation functions containing the dual operators (|2.8|) are obtained by analytic 
continuation to negative \p\ of the correlations of the "physical" states (|2^). 



There is a one-parameter family of conformally invariant boundary conditions char- 
acterized by different (generically complex) values of the boundary cosmological constant 
fx^ . In general, the boundary operators, intertwine between different boundary condi- 
tions. In the presence of boundary operators it is possible to impose different boundary 
conditions at different pieces of the boundary, each being characterised by its own value 
of fi^ . Therefore a boundary operator is not characterized only by its dimension and its 
position on the boundary, but also by the two boundary conditions it joins, characterized 
by the boundary cosmological constants fif and ■ Following [|10|, we will denote such a 

B B 

boundary operator as Bp^ ^'^ . 



If the theory is defined by the boundary term in (|2.1| ) , then the observables depend on 
the scale invariant ratios fx^ / ^/JI. For example, a disc correlation function with the bulk 

B B B B B B 

operators Vp^ , ■ ■ ■ Vp^ and the boundary operators B^j^ , 'Bq^ • • • , BgJ^ '^^ scales as 



7 = I (2 - 7str) - - V) - ^ E(l - ^2.J (2-11) 



2 ' ' "^"^ 2 

a b 



where we introduced the gravitational dimensions 



'5p = ^^, lsu = l-9 (2.12) 
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and F is some scaling function. 

The gravitational dimension 5 = 5p is related to flat dimension A = (eq. (|2.6|) ) 
by the KPZ formula 

A=^^^^. (2.13) 

1 - 7str 

The dimensions of the boundary KPZ states, or the open string states, are given by the 
same formulas, but with twice larger momentum p = 2q. 

This setting of the boundary Liouville problem is not the most general one. As we 
shall discuss below, the discrete models of quantum gravity can exhibit a critical behavior 
with anomalous scaling of the boundary cosmological constant. In this case the boundary 
interaction is e^'^. 

2. 3. Duality property of the Liouville theory 

Introduce the parametrization of the bulk and boundary cosmological constants 
through the positive constant M and the (in general complex) variable r as 

.^i^^M^ (2.14) 

TT r(i^) 

and 

^g(r) = Mcoshr. (2.15) 

TT 

Then the duality states (see fl^, eqs. (1.19) and (2.26)) that the observables are invariant 
under p, M, t ^ g,p, M, f with 

g = l/g, p = p/g, M = M\ f = gr. (2.16) 

We will assume this parametrization of n and /Us in the rest of the text. 

The meaning of the parameter r is the following. The partition function and all 
correlation functions on the disc are meromorphic functions of /x^ with a cut along the 
semi-infinite interval — oo < < —^^{0). The singularity associated with the branch 
point is a square-root one and describes the critical behavior of 2D quantum gravity with 
c = —2, in which the boundary is critical, but the world sheet has finite area. The map 
fx^ T represents the uniformization map that resolves the branch point. The two sides of 
the cut are parametrized by r = t±in, < t < oo. Note that the the boundary cosmological 
constant is an even function of the parameter r, which means that all correlation functions 
should have the symmetry r ^ — r. 
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2.4- Degenerate boundary KPZ operators 

The degenerate bulk KPZ states, or the closed string states, are labeled by the target- 
space momenta 

Prs=gr-s, r, s = l,2, ... (2.17) 

They satisfy the same fusion rule as the degenerate matter fields in the conformal field 
theory 

(r, s) * {r's') = {{r + r' -l-2m,s + s' -1- 2n)}^,^>o. 

The degenerate boundary KPZ states, or the open string states, are labeled by the 
target-space momenta 

qr — s , . 

Qrs = (2.18) 

The spacing of the degenerate target-space momenta of the open string states is thus twice 
less than that of the closed string states. 

2. 5. Functional equation for the two-point boundary correlator 

The central subject of our discussion will be the boundary two-point correlator, consid- 
ered as a function of the complex variables ri and T2 defined in ( p.l5| ) and the target-space 
momentum q 

^BMfM?BA^?A*f ^ = Dg{n,T2). (2.19) 

It satisfies the evident symmetries 

Dq{Tl,T2) = Dq{T2,Tl) = Dq{-TUT2) = D-q{Tl,T2). (2.20) 

This correlation function coincides with the coordinate-independent part of the two- 
point boundary correlator in the Liouville theory (eq. (3.18) of [1T0[| ) 

D.iru T2) = di/3\su ^2), /3 = Q/2- b\ql (2.21) 

where b = g~^^'^, Q = b + 1/b and s = r/nb (see Appendix A). 

It is also convenient to introduce the two-point correlator of the "wrongly dressed" 
states (|]8D 

Dg{ri,T2)=d{(3\si,S2), f3 = Q/2 + b\q\. (2.22) 

The "unitarity" condition satisfied by the Liouville two-point function implies a relation 
between the correlator of two "physical" KPZ states to the correlator of the "wrongly 
dressed" states (see Appendix A) 

Dq{TuT2)Dg{TuT2) = 1, (2.23) 
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where -Dg(ri,r2) is given by the analytic continuation of il>_g(ri,r2) from positive values 
of —q. 

The authors of [0 showed that the boundary two-point correlator, considered as a 
function of the complex variables ti and T2 defined in ( p.l5|) , the target-space momentum q 
and the dimensionfuU constant M defined in ( |2.14| ) , satisfies a remarkable pair of functional 
equations! 

M 

Dq{Ti + nv/g,T2) - Dq{Ti-nv/g,T2) = i-^-r~\ sinhri L) 1 (ri, r2), (2.24) 

Dq{Ti+ni,T2) - Dq{Ti-m,T2) = i—r-r sinh gn Dq_ig{Ti,r2) (2.25) 
where the multiplicative factors are given by 



These two equations are related by the duality transformation (|2.16| ). 

The equations (|2.24| )-( p.25|) hold only if all the momenta are positive. When the 
moment on the right-hand side becomes negative, eqs. ( |2.24| ) and ( p.25| ) make sense as 
analytic continuation from positive values oi q — ^g or q— ^. Using the "unitarity" relation 
( p.23|) , we obtain another pair of functional equations 

M sinhri 

Dq{Ti + iTx/g, T2) - Dqiji - iTx/g, T2) = i^Y^ 7^ 7 r, (2.27) 

Cg{q) Di_q{ri,r2) 

^ / N ^ / N sinh QfTi , 
Dqin + m, T2) - Dqin - in, T2) = i—- ^.28 



The solution of these equations is unique and reads explicitly [jTI 

Dqin, rs) = dgi\q\) Z)|,|(ri + T2)%\{n - Ts), (2.29) 
where Dqir) is given by the integral 



sinh{27Tqy/g) -^^^ 2q 



sinh(7r?/) sinh(7r?//(7) Tcy 



(2.30) 



and the normalization factor dg{q) can be found in Appendix A. 

In the rest of this paper we will study the boundary two-point correlation function 
using the the discrete formulation of quantum gravity. We will calculate explicitly the 
boundary correlators for for 

iL + l)g-l , {L + l)g-l 
q = and q = , L = l,2, ... 

and will show that they satisfy (with appropriate normalization) either (|2.24|) or ( |2.28|) . 
In both cases the answer is given by the same formula (|2.29| )- (|2.3C1| ). 



We corrected a misprint in the final formula (4.7) of |[To| ]. 
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3. Loop gas formulation of 2D quantum gravity 

3.1. The loop gas model 

The 2D quantum gravity can be constructed microscopically as strings with discrete 
target spaces. The target spaces of the rational CFT coupled to gravity are Dynkin graphs 
of ADE type These simple string theories represent ADE statistical models |[14|| 



defined on random graphs and can be realized as models of coupled random matrices [|T^ . 



A non-unitary theory with continuous spectrum of central charge c < 1 is the 0{n) model 



on a random lattice, which is realized as a one- matrix model IjT^]. All these theories are 
described in terms of nonintersecting loops of the world sheet. 

The simplest statistical model allowing interpretation in terms of the gas of loops is 
the 0{n) model coupled to 2D gravity [1^]. In this model, the matter field is defined 
as a a system of self- and mutually avoiding loops on a random planar trivalent graph. 
Each loop is taken with a length factor ne""^^"^, where (. is the length of the loop (= the 
number of links visited by the loop). The fugacity of the loops n, whose original meaning 
is the number of fiavors of the 0{n) vector field, is in fact any real number in the interval 
— 2 < n < 2. The continuous parameter n is related to the Coulomb gas coupling constant 
9 by 

n = — 2 cos Tig. (3-1) 

We are concerned only by planar trivalent graphs Q with the topology of a disc. We 
associate with such a graph the area Ag and the boundary length Lg. The area is defined 
as the number of vertices and the boundary lenght as the number of external lines of the 
graph. Then the loop gas partition function is 

2, = Q-^J■()Ag-^l^ Lg ^#ioops g-MoLioopa- ^3 2) 

graphs Q loops on Q 

where by Lioops we denoted the total length of the loops (= the number of links occupied 
by loops). 

The critical phases of the the loop gas on a random graph are qualitatively the 



same as these on a fiat hexagonal lattice [0. The dense phase is described by the vicinity 
of the critical point of Mq. In this phase the world sheet is densely covered by loops, whose 
characteristic length diverges at the critical point. The critical phase is obtained when 
both Mq and fiQ become critical. In this phase the length of the loops diverges, but the 
area that is not covered by loops also diverges. This is why the critical phase is also called 
dilute phase of the loop gas. Further, one can construct multicritical phases by adding 
more coupling constants associated with the planar graph |T^. As the loop gas model 



on a random surface can be considered as a n 7^ deformation of the ensemble of empty 
(without loops) planar graphs, its multi-critical regimes are related to the multi-critical 
regimes of the ensemble of planar graphs . 
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The critical regimes of the loop gas on a random graph are described by conformal 
field theories of matter fields with central charge c = c-^ as given by eq. (^]^), where the 
branch of the inverse function g{n) is determined by the way the criticality is achieved. 
In the dense phase < < 1, in the dilute phase 1 < g < 2, and in the /c-critical phase 
k + l<g<k + 2. 

Each value of the central charge between — oo and 1 has two realizations: one in the 
dense phase {g < 1) and the other in the dilute (or a multicritical) phase {g > 1). These 
two realizations are not equivalent, but related by duality. For example the pure gravity 
((7 = 3/2 or n = 0) is related by duality to the critical percolation problem on a random 
graph (gf = 2/3 or dense loops with n = 1). The exact solution out of criticality allows to 
determine the string susceptibility in the different phases 

(1-g, g>l; 
7str = I 1 _ i 1 < ^ < 1. (3-3) 

This means that the partition function on a surface with genus h scales as 
which in the case of the disc (genus 1/2) gives This scaling behavior matches 

with the scaling law ( p.lO|) only if g > 1. From here we conclude that the dense phase is 
described by the world-sheet action ( |2.1| ) with Liouville interaction jle^^^ instead of iie^^ ^ 
which is related to the original one by duality transformation. 

3.2. Critical behavior in presence of boundary 

In the case of a world sheet with boundary, one has to fix the boundary conditions 
both for the matter and the gravity. The boundary condition for the gravitational field 
is defined by the boundary cosmological constant, which we will denote by the letter z 
instead of /U^. This notation is inherited from the matrix model description of the loop 



gas model |T^, where the boundary cosmological constant has the meaning of a spectral 
parameter for the random matrix. 

The boundary condition for the matter field is defined by the behavior imposed on 
the vacuum loops near the boundary. The most natural boundary condition for loop gas 
model is obtained by requiring that the vacuum loops avoid the boundary. We call this 
boundary condition fixed or Dirichlet type., because in the SOS model interpretation, in 
which the loops are the domain walls separating two neighboring heights, it means that all 
points at the boundary have the same height. The partition function with this boundary 
condition satisfies a simple loop equation [jl^ . 



Another boundary condition can be prepared by allowing not only closed loops, but 
also open lines ending at the points of the boundary. The free or Neumann type boundary 
condition is defined by requiring that all links along the boundary are endpoints of such 
lines (Fig.l). This boundary condition is studied in 0. From the point of view of the 
SOS model it means that there the SOS field changes freely along the boundary. In 



8 



the continuum limit and after the Coulomb gas mapping, these two boundary conditions 
become the usual Dirichlet and Neumann boundary conditions for the gaussian field x- 




The explicit solution of the loop equations 113,3 shows that the boundary cosmological 
constant z ~ scales as /i^/^ only in the case of Dirichlet boundary in the dilute phase 
or a Neumann boundary in the dense phase. In the two other cases it scales as /U^s and 
corresponds to taking the "wrong" branch in the KPZ dressing of the identity operator 
on the boundaryil. Therefore we can identify the bulk and boundary interactions in these 
four cases as follows 



Bulk Dirichlet boundary Neumann boundary 
g>l: ^le^f n^e'^ fi^e^'^ 



Thus the duality transformation ( |2.16| ) also exchanges the Dirichlet and Neumann 
boundary conditions for the matter field, which is natural from the point of view of the 
Coulomb gas mapping. The fact that the boundary can have anomalous dimension is not 
very surprising. Take for example the topological point (7 = 1/2 in which n = and the 
world sheet has zero area. In this case the cosmological constant measures the length of 
the boundary and thus ~ /U^. 

^ This means that the Seiberg rule for choosing these branches does not necessarily hold for 
boundary operators. 
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3. 3. Star polymers in the bulk and at the boundary 



The order parameters (the magnetic operators) in the 0{n) model have a simple 
description in terms of the loop gas. The m-th magnetic operator Sl is represented as the 
source of L nonintersecting lines meeting at a point p0| . |2ll] . Following Duplantier , we 
call these operators star polymers (Fig. 2a). The correlation function of two such operators 
can be evaluated as the partition function of a network consisting of L nonintersecting lines 
tied at their extremities, moving in the sea of vacuum loops of the 0{n) model. 



In the Coulomb gas picture [|T8[ the star polymers describe point-like topological 
defects. The flat conformal dimensions of the corresponding local operators are |]2l[ 



bulk 



16 



4^7 



(3.4) 



The smallest dimension Ai is positive only in the interval —1/2 < g < 2. Outside this 
interval the propagator of the nonintersecting random walk grows with the distance be- 
tween its two extremities. The fact that two points are connected with a line leads to an 
effective repulsion between them. Such a phenomenon is typical for non-unitary theories. 
By the KPZ formula (|2.13| ), the gravitational scaling dimensions of the star operators 

are 



rbulk 



{Lg/{2iy)+jstr)/2 = ^ 



Lg g-l 



(3.5) 



These values have been confirmed by direct calculations in the matrix model [|22 





a b 
Fig. 2: Star polymers in the bulk (a) and at the boundary (b). 



The correlation function of two boundary star polymers with L lines is described 
geometrically as the partition function of the loop gas in presence of L non-intersecting 
lines whose ends meet at two boundary points (Fig. 2b). 
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The operator identification of tlie star polymers is different in tlie bulk and on the 



boundary, which is a general phenomenon in the boundary CFT [^. In the loop gas 



coupled to gravity, the boundary operators Sl have gravitational dimensions 

^bound ^L^= 5l+i,i (3.6) 

(see, for example, sect. 5 of 0.) Therefore, when considered as a boundary operators, the 
star operators Sn can be identified with the degenerate boundary fields Bg(-£,_|_i)_i. The 
fiat dimensions of these fields are 

bound L f L + 2 



Ar^'^ = 2 ( - 1 ) = ^L+1,1. (3.7) 



4. Loop equation for the disc amplitude with Dirichlet-type boundary 

In the next two sections we will derive loop equations for the boundary correlation 
functions of two star operators separating two Dirichlet-Dirichlet or Dirichlet-Neumann 
boundary conditions. These loop equations include as an ingredient the disc amplitude 
with Dirichlet-type boundary W{z). For the case of self-consistency we will derive first 
the loop equation for this amplitude. 

Let W{i) be the partition function of the loop gas on a random surface with the 
topology of a disc with fixed boundary length £. We define the disc partition function $(z) 
with boundary cosmological constant z as 

/■°° df 

^z) = / ^ e-'^ Wii). (4.1) 
Jo ^ 

Then W{_z) = ^ is the partition function with one marked point at the boundary 



oo 



W{z) = / di e"^* W{e). (4.2) 







The disc amplitude with two marked points on the boundary is given by the derivative 

dW{z) 



dz 



ide e-^^ W{e) (4.3) 







and so on. One can consider the the amplitude with k marked points and k different 
boundary cosmological constants between them. For example, the disc amplitudes with 
two marked points and two different boundary cosmological constants, zi and Z2 between 
them, is given by 

Wizu z,) = r di e--^— Wih + £,) = ^(^i)-^(^2) _ 

Jo Zi — Z2 
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Considered as a function of the complex variable Wi^z) is analytic in the complex 
plane cut along the interval — oo < z < — M along the negative real axis. The singularity 
at z = —M corns from the dominance of surfaces with divergent boundary length. The 
constant M = M{ij) is a function of the cosmological constant /U and can be calculated 
using the microscopic definition of the model It has the meaning of the effective 
boundary interaction induced by the fluctuations in the bulk. 

It is convenient to introduce, as in ( |2.15| ), the uniformization parameter r, which 
resolves the branch point singularity 

2 = McoshT. (4.5) 

Then all disc amplitudes are even analytic functions of r. 

The disc partition function W{z) satisfies a quadratic loop equation (see, for example 
§, sect 3.4) 

/d7' 1 
TT- 7 \W{z'f - 2 COS Txg W{z')W{-z')\ (4.6) 
2txi z — z 

where V{z) is a polynomial (the matrix model potential) tuned such that the ensemble of 
empty planar graphs is near the m-critical point. The contour of integration circles the 
cut in the 2;-plane along the interval — oo < z < —M on the negative half of the real axis. 
The partition function of the empty planar graphs satisfies the same equation, but without 
the last term. Therefore by subtracting this trivial piece, we can get rid of the potential 
V{z). Then equation (|4.6|) implies the following boundary condition along the real axis 



lmW{z) [W{z + iO) + W{z - iO) - 2 cos ng W{-z)] = 0. (4.7) 



In the parametrization (|4.5| ), the upper and the lower sides of the cut correspond to 
the lines t + in and r — ztt, < r < oo. Therefore ( [4.7|) can be written as the following 
finite- difference equation for the analytic function W{t) = W{z{t)) 



Its solution is 



or in terms of z 



W{t + irr) + W{t - in) - 2 cos ng W{t) = 0. (4.8) 

W{t) = cosh(£?r) (4.9) 

2 sin Tig 

W{z) = _ (^ + v^^^^)^ + (^ + V^^^)-^ ^ M 10) 

4 sin Txg ' 

where we chose the normalization for later convenience. The solution depends on the 
cosmological constant through the parameter M. If the potential V{z) is tuned near the 
m-th multicritical point of the ensemble of planar graphs, then the Coulomb coupling 
constant g should be taken in the interval m — 1 < g < m. 
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5. Boundary correlation functions of star operators 

From the operator identification ( |3.7| ) , we expect tliat tliis partition function coincides 
with the two-point function ( |2.29| ) of two boundary KPZ states with target space momenta 
q = \{{L + l)g-l) 

Wl(ti,T2) ^ D {l + i)s-i {ti,T2). (5.1) 

2 

Let us check that this is indeed the case. 

As we already pointed out, the correlation function Wl{ti, T2) of two boundary star 
operators with L lines is described geometrically as the partition function of the loop gas 
in presence of L non-intersecting lines whose ends meet at two boundary points (Fig. 2b). 
The inverse Laplace transform of Wl, which is the partition function for fixed lengths £1 
and £2 of the two boundaries 

/"OO 

WL{Ti,r2)= / rf£irf£2e-^^^i-^^^WL(£i,£2) (5.2) 
Jo 

is obtained from the disc partition function W{£) as follows. Consider the sum over all 
loop configurations with fixed lengths £'1, £']^ of the branches of the star polymers, which 
is given by the product of L + 1 disc partition functions, and then integrate with respect 
to the lengths 

/"OO 

Wl{£i,£2) = / d£[...d£'LW{£i+£[)W{£[+£'2)...W{£'L + £2). (5.3) 
^0 




Fig. 3: The loop equation for the boundary correlation function of two star polymers. 

The integral representation ( p.3| ) is equivalent to the recurrence relation between Wl 
and Wl-1 depicted in Fig. 3 

/"OO 

Wl{£iA2)= / d£W{£i+£)WL-i{£A2), Wo{£iA2) = W{£^+ £2) (5.4) 
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or, after a Laplace transformation, 

w ^ \ I W{zi) - W{z) 

Wl{zi,Z2) = <P tt- — — — Wl-i{-z,Z2) (5.5) 

J 2m z\ — z 

where the contour of integration circles the cut —oo<z<—M. This intergal equation is 
equivalent to the following condition on the discontinuity of Wl (-z) along the cut 

\mWL{-z,Z2)=\mW{-z)-WL-\{z,Z2), z>M. (5.6) 

In terms of the uniformization variable r this gives the finite difference equation 

WL{r^ + ZTT, T2) - Wl{ti - in, T2) = [W{t^ + in) - W{ti - in)]WL-i{Tu rs) (5.7) 

or, takig into account the explicit expression for W{t) 

Wl{ti + in, T2) - Wl{ti - in, T2) = ~iM^ sinh^fn Wl-i{ti, T2). (5.8) 

This equation coincides (up to a normalization factor) with the functional equation ( p.25| ) 
obtained in |T^] if we take q = ^ + constant. The constant is fixed by identifying 



Wl=o{zi,Z2)^ — — (5.9) 

Zl - Z2 

with the correlation function Dg-i{zi, Z2) of two identity operators ( [4.4|) . This gives 

^^(L+^_ L=l,2,3.... (5.10) 



which is in accord with the identification (E3.6| 
The explicit form of the solution is 



with 



A^g(L+l)-l ^ 

WL{ri,T2)= WL{Tl+T2)WLiTl-T2) (5.11) 

2 cos ng 



Lg 



m 2 



sinh \ + inkf] sinh (g^ + in ) ... sinh [g^ — 



6. The intertwiner between Dirichlet and Neumann boundary conditions 

Now let us consider another, less trivial, example of a boundary correlator, which have 
been calculated in 0. This is the correlation function 0(z,5) of two boundary changing 
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operators that intertwine between the Dirichlet and Neumann boundary conditions (see 
Fig. 10 of [|). We denote by 

2; = Mcoshr, z = M^cosh.gf (6.1) 

the cosmological constants along the Dirichlet and Neumann type boundaries, correspond- 
ingly ( see Fig. 4). 




Fig. 4- Disc amplitude with Dirichlet/ Neumann boundary condition. 



Consider first the partition function 0(£, z) with fixed length £ of the Dirichlet bound- 
ary. It satisfies the following integral equation depicted in Fig. 5 

^ /.OO 

0(£, 5) = 1 + V / W{t + 4 + ... + ~z) n{iiz)...n{inz), (6.2) 

n=o 




Fig. 5: Loop equation for the amplitude with Dirichlet/ Neumann boundary conditions 



or in terms of the Laplace transform 0(z, z) = 1 — die ^^0(£, z) 

r.r N f dz' 1 W(Z') , , 

n(z,z) = -(p , 6.3 

^ ' ^ J 2niz-z' n{-z',z) ^ ' 
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Taking the imaginary part along the cut, we get the functional equation 



ImW(z) 



or, in terms of r. 



= = (6.6) 



0(r + Z7r,r) -Q(r-Z7r,r) = ^ ^. ^, ^ ^. (6.5) 

S2(r, T) 

This equation coincides, up to a normalization factor, with the functional equation ( |2.28| ) 
for the special value q = g/4 of the target space momentum. Its solution, found in Q, is 
given by the integral ( ^.3U| ) with q = g/4. 

The Dirichlet/Neumann intertwiner Bg/4 is nothing but the gravitationally dressed 
twist operator, which should be inserted because the curvatures it/2 associated with the 
points where the Dirichlet and Neumann boundaries join. The flat dimension of this 
operator is 

(3^7 - 2)(2^-l) ^ A21A01 

For c = 1, this is the well known A = 1/16. Note also that the dimension of the Dirich- 
let/Neumann intertwiner vanishes in the case of percolation g = 2/3 and trees g = 2. 

Note that in the dual theory [g = 1/ g),m. which the Dirichlet and Neumann boundary 
conditions exchange their places, this operator should be identified as B1/4. This is con- 
sistent, since the Dirichlet-type and Neumann-type boundaries have different dimensions, 

,,B ,,1/2 ,,B ,,1/25 

^Dirichlet ~ r 5 A^Neumann ~ r 

Now it is easy to apply the method of the previous section to calculate a class of 
operators obtained by fusing the twist operator with a star polymer. Let us denote hy VLl 
the partition function of the loop gas on the disc, with two L-star polymers at the boundary, 
with Dirichlet type boundary condition on one side and Neumann type boundary condition 
on the other side. We obtain the same recurrence equation 



VIl{t + ZTT, r) - Ol(t - ^7^,r) = zsinhgfri 0L_i(ri,r2), (6.7) 

but with different initial condition Oo(r, f) = f^(r, f). This reproduces the FZZ functional 
equation for the values of the target-space momentum 

^(L+ i) - 1 
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7. Boundary correlator for generic boundary operators. 

Although obtained in quite different way, equations ( |5.8D and ( |b.5D reduce to the func- 



tional equation of [|T0| for particular values of the momentum q. These values correspond 
to the boundary changing operators that are simple to construct in the loop gas model, 
and it is plausible that the functional equation holds for any value of the charge q. 

A heuristic argument in favor of that is the following. As we saw in the previous 
section, the operator with momentum 521 = \ can be inserted by adding a line starting 
at some point of the boundary. Thus the operator boundary KPZ operator with 
momentum q can be considered as the result of the fusion of the boundary KPZ operators 
with momenta 9 — \ and q — \g- In the 2D gravity the fusion rules of KPZ operators are 
determined by the the fusion rules of their matter components. The fusion rules for the 
Liouville components follow from the requirement that the result of the fusion is again a 
marginal operator. Proceeding as in the previous section, we write the integral equation 



= i |i m^ii^^w,_,,(-.,..) (7.1) 



which leads to the functional equation 

Vrg(ri+z7r,r2)-Vrg(ri-z7r,r2) = zM3sinh^ri W ^^{rx.T^). (7.2) 



Up to a normalization factor!, this is the functional equation ( p.24|) obtained from the 



Liouville theory. The above "derivation" however applies only for operators that can be 
constructed in the loop gas model. A rigorous derivation of the functional equation in the 



general case will be presented elsewhere [24 . 



8. Concluding remarks 

We derived a functional recurrence equation for the two-point function of two discrete 
series of boundary operators in 2D quantum gravity. The recurrence equation is analogous 
to the fusion procedure with the lowest degenerate operator in the matter CFT. It is 
identical to the functional equation for the Liouville boundary two-point function obtained 
in [|10[ . We used a geometrical derivation (world-sheet surgery) of the recurrence equation, 
but it is also possible to write it as Ward identity in the corresponding matrix model. 

The method explained in this paper can be applied to arbitrary multi-point amplitudes 
and thus give an alternative procedure to derive the boundary correlation functions in 
Liouville theory [p4 |. 



^ This is the boundary "leg factor", which reflects the diflFerent normahzation of the wave 
functions in the world sheet and target space approach. 
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Since the boundary correlators in a CFT coupled to quantum gravity depend on the 
same number of variables as the boundary correlators in the corresponding CFT on the 
semi-plane, one is tempted to try to generalize the KPZ formula (p.l3| ) relating the flat 
and gravitational scaling dimensions to a one-to one correspondence between the boundary 
correlation functions on flat and random surfaces. 
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Appendix A. The boundary two-point function in Liouville theory 

The authors of [|TD[ considered pure boundary Liouville theory with action 



= y + '"^'''^) + e'-^ (A.l) 

where the background charge is introduced by the asymptotics 

(Piz,z) ^ -Qlog{zz), Q = b + l/b. (A.2) 
The Liouville bulk and boundary vertex operators are defined as 

Considered as a function of the boundary cosmological constant fx^ , the disc partition 
function is analytic in the complex plane cut along the interval 



—00 < < — a//U sinyrfe^. 
The branch point singularity is resolved by the uniformization map fx^ s defined as 



jj,^ = \/ \i sin 7r62 cosh(7r6s). (A. 4) 

The boundary 2-point function 

'Bf^"{x)Bf^"{^)) = rf(/3|^f l^r^^'^ (A.5) 

18 



is given by eq. (3.18) in [10|. It is related to the correlation function of two boundary KPZ 
states in 2D quantum gravity as 



In the notations used in this paper 

i^ = h(l 



B _ j Dq{Ti,T2) if /3 < f 



9 = ^^ 



P 



M 



Q-2a 



a//x sin7r6^ 

TV—— -TT^ 



r = nbs, 

Q-2(3 
2b 



r(i-62) 

the coordinate-independent part of the boundary two-point function reads 

d{f3\iSl,S2) = dg{q) Qq{Ti +T2)Gq{Tl - T2) , 

where 



log^g(r) 



1 fdy 

27 y 

—00 



sinh{27rqy/g) 
sinh(7ry) sinh^Tcy / g) 



2q 

Txy 



and the function G is defined by (6 = 1/ y^) 

G{x + b) = ^^^r(6x)G(x), G(x + 1/6) = ^^-^^^r{x/b)G{x). 



/2tt 



27r 

The function Qq{r) satisfies the unitarity condition 

^,(r)^_,(r) = 1. 

and the recurrence equations 



Qg+igir) = cosh [ — ±inq)gq{TT ^tt 



1 



(A.6) 



(A.7) 



(A.8) 
(A.9) 

(A.IO) 



(A.ll) 



(A.12) 

(A.13) 
(A.14) 

Qq+^ir) = cosh ^- ± ZTT^ J Gqir T ^^7r) (A.15) 

which determine it up to a multiplicative constant. For boundary operators corresponding 
to degenerate KPZ states q = \{rg — s) the function Qq takes a simple form 



Q rg-s (r ) 



j-jr-l ^^^^ gr+^^[^r-2k-l)g-s] 

n =; cosh / — ^ 



(A.16) 
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